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1 Introduction 



Saturation effects in the parton distributions for a nucleon, that could finally 
lead to unitarization of the cross sections, have been thoroughly discussed 
in the past. In the small Bjorken x region these effects can be accounted 
for by the introduction of non-linear terms in the evolution equation for the 
gluon density p. The study of non-linear evolution equations began twenty 
years ago and gave rise, in the following years, to different approaches to the 
small-a; region [3 El • More recently, generalizations of non-linear evolution 
equations have been proposed with different physical motivations [H El E] 
and screening effects have been incorporated in the theoretical framework in 
different ways HllHlli. 

The approaches quoted above find interesting applications in the interac- 
tion of ultra- high energy neutrinos with nucleons and nuclei fOl ^1 ^1 ^] ■ 
While there is no general consensus on the importance of screening effects, 
one expects, at any rate, that the linear QCD evolution of parton distribu- 
tion functions will be tamed in the very small-x region. Predictions for the 
cross section of the neutrino- nucleon interaction based on DGLAP jHl ^] or 
BFKL jinj equations show a power increase with energy ^3 ^1 ^\ that will 
finally violate the Froissart bound. We notice that the problem these papers 
address is rather intricate, since it requires a complete knowledge of nucleon 
structure functions in both variables x and Q^. At ultra- high energies the 
integrals, giving the cross section in terms of structure functions, cover the 
whole permissible range of values for these variables and explore extreme re- 
gions of the (Q^, x) phase space, where non-accelerator data exist. Moreover, 
the mathematical complexity inherent in the solution of non-linear evolution 
equations may conceal the physical essence of the problem. 

Simplifications are well possible if we limit ourselves to the small-x region 
only, with a warning about the consequences on the value of the integrals 
giving the total cross section. For example, a power-like behavior in x of the 
parton distributions is a simple solution of the DGLAP dynamics at next- 
to- leading order on the region x < 0.1 and for large values, where the 
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hard initial condition x~ ^ constant can be applied for parton distribution 
functions ^7j . Due to the asymptotic behavior of the DGLAP evolution ^3] , 
small-x data can also be interpreted in terms of the "double asymptotic 
scaling" [201 1^ thai provides an explicit solution to the problem. The 
straightforward application of these approximation schemes to the non-linear 
evolution is however questionable. Moreover, it would be rather difficult to 
estimate the error induced on the cross section. 

In a series of papers [2S1 121j a new set of evolution equations was sug- 
gested that includes parton recombination. These new equations are derived 
in the leading logarithmic (Q^) approximation and differ from the traditional 
ones mi2], that rely on the double logarithmic approximation (DLL A). DLLA 
means that only those terms in the splitting functions that generate large log- 
arithms in X are important. In other words, in the DLLA a diagram consists 
of gluon ladders and any transition from gluon to quark is suppressed in this 
approximation. It becomes difficult however to reconcile small-x approxi- 
mations with the twist-4 gluon recombination corrections of Ref. |24|. As 
emphasized in Ref. the twist-4 coefficient function, driven by a two- 
particle gluon distribution, cannot be simplified by using DLLA or other 
small-x approximations. 

There is, however, the possibility to replace, at small x, the convolu- 
tion of two functions by a simple product. The method, introduced in 
Refs. |2ni nil I2H1 121]; allows a correct treatment of the non-singular part 
of parton distributions and has found numerous applications [23 IHOl 1^ - 
This approach, to be described at length in the following, has been applied 
also to the evaluation of the contributions from higher-twist operators of the 
Wilson operator product expansion jHIj. The accuracy of this method, when 
applied to the modified DGLAP equations of Ref. 124,) can be verified a 
posteriori with a suitable computer code. 

In this paper we study the contributions of the twist-4 gluon recombi- 
nation corrections to a previous twist-2 calculation [30]. In Ref. we 
estimated the ultra-high energy neutrino-nucleon cross section in the ap- 
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proach of Ref. The corrections introduced by the presence of non-hnear 
terms in the evolution equation, as given in Ref. [23) will change sensibly 
our previous estimate. 

In the next Section we will discuss the non-perturbative input and present 
an analytical form for the evolution. Some restrictions on the value of the 
parameters of the input distributions, present in Ref. jSU], are relaxed in this 
new formulation. However, gluon recombination corrections will be estimated 
in the particular case of starting flat initial conditions. In Sections El results 
are presented for the structure function FJ^^i^x^Q'^) and for the charged- 
current neutrino-nucleon total cross section. The asymptotic behavior of 
this cross section will be also discussed. In the Appendices the proof of the 
relevant analytical results is presented. 

2 evolution 

As in our previous paper [30] , we choose a soft non-perturbative input based 
on analyses of the nucleon structure functions [S21 IHHl IHH IHEl IHEl IHZI • If we 

denote by Q^) the sea quark distribution xS'(x, Q^) and by fg{x,Q'^) 
the gluon distribution xG{x, Q"^), that is if we put 

Q') ^ xSix, g2), Q') ^ xGix, Q') , (1) 

our soft non-perturbative input can be written in the form 

fa{x, Qq) = 

where Aa, Ba and i^iQl) are unknown parameters to be determined from 
data. Throughout this calculation at small x we will ignore the non-singlet 
quark component and limit ourselves to the leading order (LO) of perturba- 
tion theory. 

The factor ( 1 — has been treated in detail in Ref. [301 ^^"^ we neglect 

it in the following. At the end we will take it into account by multiplying the 
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resulting structure function, F2^{x, Q"^) = fg{x, Q^) in the case of neutrino- 
nucleon DIS, by an effective large-x behavior 
Furthermore, as in Ref. jSUj, we define 
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with constant u. 
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and the Ball-Forte scaling variables 



a = 2y —dggt\n{l/x) 



P 



-dggt 



a 



(3) 



where dgg = —12/j3o and /3o = H 



ln(l/x) 21n(l/x) ' 

2//3, with / the number of flavors, is 
the LO coefficient of the QCD /5-function (in units of — IGvr^). For brevity, 
we introduce the notation = 1 + 20//(27/?o), d^{l) = 16//(27/?o)- At 

the LO, a,{Q') = 47r/{(3, In(QVAio))- 



2.1 DGLAP evolution 

In Ref. ISO] the results of Ref. jSHj were used in order to obtain an approximate 
evolution, in LO perturbation theory, under the conditions Aa ^ Ba (here 
and in the following the index a stands for q ot g), so that no interference 
appears in the evolution of the coefficients multiplying the different powers 
of the logarithm. Since here we relax these conditions, it is interesting to 
present the new expressions for fg and fg. 

The method of solution adopted in Ref. [201 can be summarized as fol- 
lows. Starting from the exact solution in the moment space, the anomalous 
dimensions and the coefficient functions are expanded in the neighborhood 
of n = 1. The singular part, when — > 1, leads to Bessel functions but, in 
order to achieve the accuracy 0(p), also the regular part must be properly 
taken into account in the inverse Mellin transform. 

By analogy with Ref. it is possible to obtain the small-x asymptotic 
results for parton distribution functions (PDFs) and the F2 structure function 
at LO of perturbation theory by setting 

fa{x,Q') = f^{x,Q') + f-{x,Q'), (4) 
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where 



X 



. p-d-Wt 



+ 0{x) 



A+p/i(a) + i?+/o(a)l -e-^+W*- fl + 0(p) 



(5) 

(6) 
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-d+{i)t 



1 + 0(p) . (7) 



Here, In{z) are modified Bessel functions and tlie coefficients A"^ and are 
defined as 
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5 9 V 9 ^ ' ' 



i^B. + ^Bg, , 
4 ^ 2 

A„ + -Ag + 

9 * 27 



9 V ^ 9 ^ 6 



A„ + -A, - - ( 1 + 1^ 1 5„ - 



27^-9 3 ' ' 
27/ ^ 243 ^ 



(8) 



2.2 Shadowing and anti-shad owing corrections 

According to Ref. [24j, sea quark and gluon distributions are modified by 
the introduction of gluon recombination as stated by the following modified 
DGLAP equations: 



c/lng2 
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-K2 



(9) 



Here Pab'^x) are the Altarelli-Parisi kernels, ® stands for the Melhn convo- 
lution, defined as 



In Ref. it was set Ki = K2 = K and from a fit of the HERA data the 
resulting value for K in the evolution equations of Ref. ^Ij turned out to be 
very small, K = 0.0014. We will come back to this point in Subsection 13.11 
The introduction of two parameters, Ki and K2, allows a clearer check of 
the importance of the anti-shadowing, with respect to the shadowing contri- 
bution. Moreover, they give the possibility to relate models introduced in 
Refs. 123 121 and m ig (see Subsection EH). 

Since Ki and K2 are expected to be small numbers, the solution of Eq. © 
can be written as 




and 




(10) 



(11) 



//""(x,Q2) = /.(x,g2)+T,(x,g2), 



(12) 



where 




(13) 



and 



dT,{x,Q^) 
cilnQ2 



b=q,g 



(14) 
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with 



Finally, the sea quark distribution is 



xS{x,Q'') = ff^\x,Q''). 



(15) 



2.3 Small-x solution of the complete equations 

The solution of Eq. (|T3|l with the boundary condition (j2)) has been found 
already and is expressed in Eqs. (jD)-(|7j). This result has been obtained di- 
rectly from the corresponding solution in the moment space (see Ref. |29j). 
In order to simplify the solution of the complete equations, we set in the 
following Bq = Bg = 0. Then our solution ®-(IZj) assumes the simple form 



A>Ji(a)-e-^+«*-(l + 0{p) 



f:{x,Q') = A+/o(a)-e-^+«*-(l + 0(p)) . 



(16) 
(17) 

(18) 

Equation ()14|) can be rewritten in the moment space, using the Mellin trans- 
form defined as ^ 

M{n) = I dxx^-^M{x) . 
Jo 

In the leading logarithmic approximation we obtain 

dTa{n,Q^) 



dt 



J2 dabTt{n,Q^) + ra{n, Q 

b=q,9 



(19) 



where dab{n) = 7^°''(?t,)/(2/5o) is the ratio between the anomalous dimension 
7^°^ and twice /3o and ra{n, Q"^) = 47ri?a(n, Q'^)/Po, with Ra{n, Q"^) the Mellin 
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moment of Ra{x, Q^). The solution of Eq. (|19|) can be written in the form 
Un,Q') = T^{n,Q') + T~{n,Q') , 

/oo 
dw e'^^""-'^ J2 'aMn, M') , (20) 



b=q,g 



where 



w = In 



d±{n) are the eigenvalues of the dab matrix and e^j,(n) are related to the 
components of the eigenvectors of the same matrix. Explicitly, we have 



dw e'^+("'~*) 

dw 6"^+ 



(1 - r])rq{n, M^) - firg{n, M^) 
(1 - ri)rg{n, M^) + erq{n, M^) 
r]rg{n,M^) - erg{n,M'^) 



(21) 



where e^^, have been expressed in terms of rj, rj and e, which take the following 
values for n ^ 1: 



ri = 1 In — 1, 



--in -I] 
9^ 
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(22) 



The evaluation of the Mellin moments of Ra{x^Q^) requires particular 
care, since the integrals appearing in Eq. (|T5|) are not exact Mellin convolu- 
tions. If we make the position 



M{n\l/2) 



1/2 



dxx''~'^M{x) 



we have, in addition to the usual Mellin convolution 



dx 



' dy 



X 



Ml - M2(y) = Mi(n)M2(n) 



(23) 
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also the following one: 

x^'-'dx ^Mi MM = 2"-iMi(n)M2(n|l/2) , (24) 

Jo Jx/2 y \^yj 

where the definition Mi{x) = Mi{2x) is used hereafter. Now, in order to 
calculate the Mellin moments of Raix^Q"^) in Eq. (fT3jl . we should find a 
relation between the "Mellin transforms" f2{n) and f2{n\l/2), defined as 

h{n) = [' dxx--'f^{x) , (25) 
Jo 

h{n\l/2) = j'^\xx--^fl{x). (26) 

This relation is 

/2(n|l/2) = /2(n) + 0((n-l)2), (27) 

the proof being given in Appendix A. 

Then, the contributions to the Mellin transform of i?a(x, Q^) which are 
regular for n — > 1 assume the form (the upper index (r), here and in the 
following, stands for "regular contribution") 

' F^l\n) - Fil\n)\ Un) , (28) 

for the coefficient of Ki and 

Fil\n)Un), (29) 
for the coefficient of i^2, where (see also Eqs. and (j211)) 

F^in) = [\xx--'F(l\x), 
Jo 

Fil\n) = fdxx-'F(l>{x)=[ dxx--'Fil\2x), (30) 
JO Jo 

since by definition Fal\x) = Fag\x/2). 
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Note that Fqg\x) = Fqg{x), because the Melhn moments of Fqg{x) have 



no singular part, and 

27 
64 

By performing the Melhn integrals in Eqs. ()30|) . we find 

1813 



-99x - 4x(2 - x) ( 34 - 33x + 16x^ - 4x^ 



(31) 



2880 ' 
131 

180 ' 
31977 

320" ' 
23877 

160" ■ 

In this way we obtain from Eqs. p5|) . (j^Hj) and (j^^ 



Rq{n 1,M2) 



Rg{n 1,M2) 



2880 



2880 



M2 
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X 



32 
M2 



297, „ 1753\ 
In 2 — 



320 y 



Ki 



297 1 3553\ 
32 n - 1 ~ ^20" J 



(32) 



(33) 

K2 

(34) 



and, from Eqs. (pTjl 
T±(n^l,Q2) = _ 



(M2) 



2xxd±(n^l) 



(a,(M2))''±("^^)-' ■ 
(35) 



Here C^g{n — > 1) are the "coefficient functions" of the twist-4 corrections. 



because 



Cagin^l) = C+(n^l) + C-rn^l), 



(36) 



i.e. they are the coefficients in front of the Mellin moments of the function 
fg{x). They can be written as 



ag 



ag 



(37) 
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and the definitions C'^^^in 1) = C^i^ {i = 1,2) have been introduced for 
the sake of simphcity. The coefficients C*^ are given in Appendix B. 

The solution ()35p in the moment space can be transformed to the x-space. 
Note that the product of moments in Eq. (jH^j) leads to the convolution 

Mi(n) ■ M2(n) ^ M^ (-] M^iy) = M,{x) ® M^ix) (38) 

Jx y \yj 

in the x-space. As in the case of the moment space, Ta{x, Q"^) can be repre- 
sented as the combination of the "+" and "— " components: 

Ux,Q')=T^{x,Q') + T-ix,Q') , (39) 

which can be obtained from the corresponding components in Eq. (|35|) by an 
inverse Mellin transformation. 

As for the "— " component, we note that the value d^{n) does not contain 
any singularity for n — 1, hence (hereafter v = w — t) 

gd_(n). _ e'^'^^>6{l-x) (40) 

(here 5(1 — x) is the Dirac 5-function), so that 

e'^-(^)(-*) Mn 1, M') C-^ e'-'-'^^-'^ /|(x, M') . (41) 

As for the "+" component, we have d^{n) = d^/ (n— l) + (i+(n), with < 0, 
hence 



M- 



]^ d+v/{n-l) ^ ]_ {d+'^Y 

- 1 in- 

oo / 1 \ ^ 



fc=0 

J+v/(n-l) ^ ^ }_ jd+v)'' 



k\ in -If 

oo _ / 1 \ 



k=0 

S{i-x)-pMa) , (43) 
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where 



P = 
P = 



21n(l/x)' 



21n(l/a;)' 
and w — t when Q"^ M^, i.e. 



a — a with t —>■ w , 
a = a with t ^ {w — t) 



w 



w — t 



In 
In 



In 
In 



V ln(Qg/A 



2 "I 
LO) 



/ln(MVA2 



As a consequence 

C+ (n = 1) e'^+("'-*)/("-i) e^+(i)("'-*) /2(n ^ 1, M^) 



M- 



M- 



1 



C+ (n = 1)1 e'^+(— e'^+(i)(^-*) f2{n ^ 1, M 



where 



r 

y 



S{l-y)-p{y)Ma{y)) fJx/y,M') 



5{l-x)-pMa)\®f^{x,M'), 
' - Maiy)) f^{x/y, M') ^ ^a) ® f^{x, M') . 

y 



Thus, 

T-{x,Q') 



X 



d-(l)-2 ' 
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X 
X 

The above formulas complete our calculation of the gluon recombination 
terms. 

2.4 Evaluation of Fi and F2 at 0{p) 

The functions Fi and F2 can be evaluated approximately. Note that we can 
write 

F2 = F,- + 2F+- + F++ , (55) 

according to the decomposition (see Eq. (jl))) 

/|(x, M') = (/-(x, M')f + 2/;(x, M')f;{x, M') + (//(x, M^))^ . 

The details of the complete calculation, at order 0(p) and with Bq = Bg = 0, 
will be given in Appendix C. Here we present only the results. 
For the " " component we find 

Ff- = (A;)Vo(a)-e-^'^-«-(l + 0(p)), (56) 

F,- = (A;)2iji(a)-e-2'^-«-(l + 0(p)) , (57) 

while the result for the "H — " component is 

F+- ^ -pJi(a)® (/;(x,M2)/-(x,M2)) 
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Fi 



(a,(M2))^+W-2 



(53) 



Cgg^F2+Cgg ^1 



2N\a!+{l) 



(«s(M2)) 



2NNd+(l)-2 



(54) 



A;AjJo(a) ■e-(''-«+'^+«)"'(^l + 0{p) 
Ma)®(f^{x,M')f;{x,M^) 



P 



(58) 



(59) 



The term corresponding to the component "++" requires a more involved 
treatment. As shown in Appendix C, one obtains 

= 5(l-x)-pJi(a)l ® f/+(x,M2)'' 



Ft 



and 



(A:)^/o(V/^)/o(^V/^)-e-^^+«--fl + 0{p) 



Ma)0{f+ix,M') 



(60) 



= (A+)2-^y/^[./o(.V/^)/i(V/^) 

- -zIo{zy'/')hi-zy'/')] ■ e'^^-W- ■ (l + 0{p')) , (61) 

where y = — log(a;), z = a^l'^ + % with a = —d+{3w + t) and /? = 

—d^(t — w), and z is the complex conjugate of z. 

These approximations do not make the integrations in Eqs. ()53|) - ()54p 
avoidable, but greatly simplify the analytic structure of the answer. 



3 The neutrino-nucleon cross section 
3.1 The fitting procedure 

We consider the parton distribution functions from the ZEUS Collabora- 
tion [38j in the region 2.5 GeV^ < < 20 GeV^, for values of x in the range 
10~^ < X < 5 X 10"'^, where the ZEUS NLO fit favorably compares with ex- 
isting HERA data. From the parton distributions it is possible to reconstruct 
the isoscalar structure function 2xFi for neutrino-nucleon scattering 

2xF^^{x) ~ 2xF^^{x) 
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^ xu{x) + xu{x) + xd{x) + xd{x) + 2xs{x) + 2xc{x) + . . . (62) 

where . . . stands for b and t quarks and we have assumed s = s and c = c. 
At the LO -^2(2;) = 2xFi{x) and the difference between the two structure 
functions decreases when increases. Parton distribution functions are 
needed at any rate; the knowledge of the experimental must in fact be 
supplemented by PDFs since 

Fr=^-^F^'' + %s-xc) + ... 
5 5 

For the values under consideration, there are no CCFR data [HU] below 
X = 0.0125 and the comparison between our fit and F2 measurements in the 
process (z/^ + z/^) + nucleon — > {fi~ + + X is impossible. A single high 
energy HERA data point [IHl for the process e~p —>■ u^X, with p± > 25 GeV, 
gives cr = 55 ± 15 pb at ^/s = 296 GeV. This means for the neutrino-nucleon 
cross section a value of (2.0 ± 0.55) x 10"^^ cm^ at ^/s = 306.4 GeV. 

The possibility to test other saturation models, such as that of Ref . PP , 
relies on the presence of two different coefficients Ki and K2 for shadowing 
and antishadowing contributions. In the following, we set Ki = K2 = K, 
as in Ref. [21] , thus reducing the number of parameters in our approach and 
the related errors. Then, the free parameters become A^q, Qo: K. 
We choose to fix / = 4, as in Ref. [31J, and A^o = 0.19 GeV. The remaining 
parameters should be determined by a fit. There is one more parameter to 
be considered, namely the "mean power" u of the factor (1 — x) present in 
F2^{x, Q'^). The term (1 — xY will be considered, for the sake of simplicity, 
only at the end of the calculation and will not evolve in our model. This 
will surely affect the result, but, since we are interested in the total cross 
section at very high energy, and hence at (x) very small, the error should not 
influence too much the result. 

The values of the parameters Qq and K can be estimated from other 
sources. In the paper by W. Zhu et al. [21] a fit to the HERA small-x data 
for has been performed starting from GRV-like input distributions 
at Qo = 0.34 GeV^. MD-DGLAP evolution equations Plj determine new pa- 
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rameters for the sea quark distributions at Ql, with respect to GRV98LO jH], 
and give for the non-hnear coefficient K the value K = 0.0014. 

On the other hand, the fit of Ref. [HI] imposes a constraint on the param- 
eter Qq. This constraint originates from the conditions of apphcabihty of the 
"generahzed double asymptotic scaling", that gives a satisfactory descrip- 
tion of the experimental data in a region > QLt) where Qc^t is a cutoff 
larger than Qo- For small x and large values it is reasonable to neglect 
the valence quarks and the disagreement at small is a consequence of 
this approximation. The value of Ql in this approach, that we follow in our 
paper, turns out to be rather small and, for the LO fit, it is approximately 
Ql ~ 0.3 ^ 0.4 GeV^. This value of Ql has been obtained by fitting HERA 
data with higher-twist contribution evaluated in the renormalon model [32] ■ 

At this point it can be useful to remember that the purpose of this calcu- 
lation is to reproduce the neutrino-nucleon cross section at very high energy, 
where presumably valence quarks do not contribute. Looking at Ref. [TH] . 
in particular at the Figure 2 in this paper, one sees that the valence quark 
contribution to the total uN charged- current cross section is larger than the 
u and d sea quarks contribution up to the energy ~ 10^ GeV. Hence we 
cannot rely on the HI Collaboration data point at E'jy ~ 5 x lO'* GeV and we 
expect that our model can be trusted only at higher energies. 

3.2 Results 

We have reconstructed the structure function F2^{x,Q'^) from the ZEUS 
PDFs j2H] and performed a fit to the asymptotic formula for F2^{x,Q'^) = 
where has been defined in Eq. (fT^ and the ingre- 

dients for calculating it have been given in Section |2l Basing ourselves on 
the arguments of Ref. ^Tj, we impose an upper bound on Ql, Ql < 0.45 
GeV^, and consider separately two possible values for the mean power u of 
the factor {1 — x): z/ = 4 or z/ = 5. The lowest x^/d.o.f for this constrained 
fit has been obtained having chosen Ql = 0.45 GeV^ and z/ = 4, with the 
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result (xVd-o.f.=0.553) 



Aq = 1.002(39) 
Ag = 0.565(29) 

K = 0.0130(49) . (63) 

There are several reasons why the parameter K in ()63p is much larger 
than the corresponding result in Ref. First of all, the change from the 
leptonic structure function of Ref. |25 to the neutrino one requires a factor 
near to 18/5 that affects also the coefficient K. In addition, our fit shows that 
there is a strong correlation between the starting value of Q^, i.e. Ql, and 
the value of K, which turns out to increase with Qq. These reasons suffice 
to say that, within errors, our result is compatible with that of Ref. 

In order to test our solution, we have calculated F2'^{x,Q^) from the 
ZEUS PDFs |HH] at = 50 GeV^ and compared with our theoretical result. 
The percentage error at x = 10~^ is nearly 2% and is of the same order of 
magnitude in the whole range of x: 10~^ < x < 5 x 10^^ (see Figure 1). 

Another test of this model regards the slope dF2^{x,Q'^)/dlnQ'^. As 
noticed in Ref. |2S] , higher- twist effects are more easily revealed in the slope 
than in F2{x, Q^). Since in Ref. [23] it was stated that one cannot rely on the 
DLLA for the calculation of the twist-4 contributions, a critical examination 
of our approximation becomes important. It is easy to write in our approach 
(see Subsections 12.21 and 12. 3j) 

dFr{x,Q') _ df^-'Kx,Q') df,{x,Q') ^ a^^ (-^lf{x,Q')] . (64) 



c/lng2 d\nQ'^ c/lnQ2 Q2 \^ ^2'^ ^ 

In Figure 2 we show the behavior of the slope as a function of x for = 10 
GeV^ in the range 10"^'' < a; < 10"'^ according to our theoretical results, 
with the parameters given in Eqs. (j63|) . A similar behavior is observed for 
other values of the range 2.5 GeV^ < Q'^ < 20 GeV^. From this 

Figure it is possible to see the higher-twist effects and to have the idea of the 
uncertainty resulting from the the errors on the fitted parameters. 
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Finally, we consider the total cross section. With the usual notation, we 
can write for charged-current neutrino interactions 



da y _ GlME^ ( M, 



/'2 



dxdyj TT VQ^ + ^h'. 

X {^-^^^^p^F^M^) + (l - I) yxFr(x,g^)} .(65) 

For anti-neutrino charged-current processes one must change the sign in front 
of F^^{x, Q"^). In Eq. Gp is the Fermi constant, M is the nucleon mass 
and the variable y is related to the Bjorken x through the relation 

= ~ 91 (m) 

^ x{s-M^) ~ xs ■ ^ ' 

The laboratory neutrino energy = {s - M2)/(2M) is approximately 
s/{2M) in the energy region of interest. The F2 contribution to the total 
cross section can be written in the form 

-uN _ ^ + 



2 



r ,r,2( Ml y rfa; 1 + (1 _ g2/(^,))2 



(67) 

and the formulas obtained for the FJ^^^x^Q"^) allow the evaluation of o"^^ 
that, for large s, is a good approximation for the total vN charged-current 
cross section. 

In Figure 3 we show the behavior of a'^^ as a function of s, in the range 
10^ GeV^ < s < 10^^ GeV^. The relative error is 4.4% at s = 10^ GeV^ and 
14.0% at s = 10^^ GeV^. For comparison we draw on the same plot the results 
obtained in Ref. and in Ref. ^Hl as well as the HERA measurement at 
^/s = 306.42 GeV. At this value of s our estimate of the cross section is much 
lower than the HI data point. However, this is not surprising in view of the 
fact that we neglected valence quarks, which still can play a residual role at 
these relatively small values of s. For larger values of s, our results nicely 
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compare with those of Refs. fHl showing good agreement till s ~ 10^^ 
GeV^. For even larger values of s the effect of higher-twist starts to become 
visible. 

These findings support the conclusion that our approach, based on an 
asymptotic calculation, succeeded in singling out the relevant part of the 
cross section. 

In order to make the higher-twist effects more visible, we determined also 
a"^ without considering the recombination terms in the evolution equations, 
i.e. with K = from the beginning. In this case the fitted parameters 
turn to be Ag = 1.040(36) and Ag = 0.548(28) (xVd.o.f.=0.665). In Fig- 
ure 4 we compare our results for a'^'^ with and without the inclusion of the 
recombination effect. The higher-twist effects become visible for s = 10^^ 
GeV2. 

4 Conclusions 

Twist-4 corrections to the structure function F2 have been estimated at small 
X in leading order QCD following the method developed in Ref. |2ni- This 
estimate leads to an analytical parametrization for the gluon recombination 
effects and completes, in this respect, the program outlined in our previous 
paper PU], where an approximate QCD evolution at twist-2 was presented. 

Our approach to the saturation phenomenon follows the scheme pro- 
posed in Ref. j2S]- The non-linear evolution equations we study are the 
MD-DGLAP equations [21], where momentum conservation gives rise to an 
antishadowing term that influences appreciably the screening effects. A com- 
pact and analytical solution of these equations at small x is possible only if 
some conditions are satisfied. The most relevant of these conditions regards 
the input parton distributions that must be fiat. The proof of Eq. (jlTj), which 
allows a simple treatment of the gluon recombination terms, requires this as- 
sumption. Moreover, the kernels in the MD-DGLAP equations are given in 
leading order of perturbation theory and this approximation reflects on our 
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approach. In spite of these assumptions and approximations, the results are 
satisfactory. Examples of the Q"^ evolution and of the behavior of the slope 
dF2/d\n.Q'^, shown in Figures 1 and 2, make us confident on the accuracy of 
the method. 

The most interesting application of our simplified formulas is, in our opin- 
ion, the study of the interaction of ultra-high energy cosmic neutrinos with 
nucleons. This process can probe an energy region far beyond the largest 
energy reached by existing accelerators. The simplified and reasonable ex- 
pressions for the structure function F^^, disposable in this paper, renders 
the evaluation of the neutrino-nucleon cross section and the estimate of the 
twist-4 contributions simpler and more transparent. Our result for the cross 
section, shown in Figure 3, is in perfect agreement with the calculations of 
other models and, only above s = 10^^ GcV^ {E^, ~ 5.3 x 10^^ GeV), the 
effect of twist-4 gluon recombination becomes visible. The Froissart limit 
will be eventually satisfied, but in a region where the complete solution for 
the gluon distribution becomes necessary. 

It is well possible, at such high energies, to simplify further the integrals 
leading to the cross section. Work on this problem is in progress. 
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Figure 1: Comparison at = 50 GeV^ between F2'^{x, Q"^) as obtained from the 
ZEUS PDFs 38 and from our theoretical calculation (the vertical bars represent 
the uncertainties coming from the error in the fitted parameters). The values of 
the parameters of Eqs. ()63|) entering our results were obtained by a fit performed 
in the region 2.5 GeV^ < Q"^ < 20 GeV^, 10^'' < x < 5 x 10"^. 
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Figure 2: Slope dF^^ {x,Q'^)/d\nQ'^ at = iq GeV^ according to our results. 
The uncertainties on the data come from the errors in the values of the fitted 
parameters of Eqs. 
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Figure 3: Cross section a as a function of s according to our results (filled 
diamonds). Data do not include the uncertainties coming from the errors on 
the values of the fitted parameters of Eqs. (|63j) (see, instead, Figure 4). For 
comparison, data from Ref. [IHj (open circles) and from Ref. UHI (open triangles) 
are also shown. The isolated point at ^/s = 306.42 GeV (filled squares) represents 
the HERA measurement. 
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Figure 4: Cross section a" as a function of s according to our results (filled 
diamonds, same as Figure 3, but with the inclusion of error bars) and without the 
contribution from the recombination (open squares). 
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A Appendix: Proof of Eq. ( 1271 ) 

We consider separately the three terms coming from the decomposition /| = 
ifgf + 2/3"/+ + (/+)^ which will be named in the following "++", "+-" 

and " " components, respectively. Moreover, we will use the notation ^ to 

denote the Mellin transform and to denote the inverse Mellin transform, 
a) For the " " component (hereafter A is arbitrary) we get 



y ' " Jx y X 
1 , .\ , „ , M ( 1 



ln--lnAj (/-(x))^ ^ (___lnA)/2-H, (A.l) 



where we have used the definition 



n\n)= l\xx-~'fUx)mx) (^,j = ±). (A.2) 



The same argument gives 



Jx/2 y Jx/2 y X 



r/2 y Jx/ 

ln^-lnA^(/;(x)f ^ (^_i_ + ln2-lnA^/2--H.(A.3) 
Then, from Eqs. ()A.1|) and ()A.3|) we obtain 

r ^(/,-(2/))' = ln2(/-(x))2 ^ ln2frin). (A.4) 

Jx/2 y 

b) For the "H — " component (hereafter z = \n{l/y) and A = \dgg\t) we get 

Jx y Jx y J In A ,_n 



^ k\(k + r 

fc=0 ^ ' 



In- -(In A) 
X I 



k+l 



(k\ 

k=0 ^ 

P PA 



i/i(a)-(l + 0(p2)) ^ (-l^ + 0(n-l))/+-(n), (A.5) 
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where 

Pa = p\x-^i/A, (^A = (^\x^i/A ■ (A. 6) 

Similarly, we obtain 

r/2 y 



/■^/^ dv 11 

/ - /;(?/)/,- (y) - -h{cr,) - -h{aA) , (A.7) 
Jx/2 y P2 Pa 



where 
Note that 



P2 = p\x-*x/2, (^2 = (y\x^xl2 ■ (A. 8) 



A' 



= -/i(or) + ln2Jo(a) ■ (l + 0(p)) . (A.9) 



P 

Then we arrive at 



jf^{y)f;{y) - -^h{a)+\n2Io{a) -{l + Oipi 

^ (;^ + ln2 + O(n-l)^/2+-(n).(A.10) 
Finally, from Eqs. ()A.5|1 and ()A.10|) . we obtain 

r -f^iy)f;iy) = ln2/;(x)/;(a:)-(l + 0(p)) 

Jx/2 y ^ ^ 

^ (ln2 + 0(n-l)) /2+-(n). (A.ll) 
c) For the "++" component we get 



rVA , .l/A , .ln(l/x) 

/ -Utiy)?- ^iMy))= dzi^ 

Jx y Jx y JlnA 



a{z)) 



-h{a)Io{a) . (l + 0(p2)) - / dzP,{a{z)) , (A.12) 
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where we have used the integration by parts and the results from the previous 
case b). 

Since, for a ^ oo, Ii{cr) = /o(o") ■ (1 + 0{p)), we have 

f^^j^f^iy))" - ^A(^)/oM-(i + o(p^)) 

- (^^ + 0(n-l)^/++H. (A.13) 



In the same way, we obtain 

-l/A 
r/2 



[ -U:{y)f - ^/i(^2)/o(^2)-(i + o(/f 
Jx/2 y 2p2 V ^ 

= ^/i(a)/o(a) + ln2/o2(a) -(l + Olp)) 



^ (;^ + ln2 + 0(n-l)j/++H (A.14) 

and, from Eqs. (IA.13I1 and (TOil) . 

= ln2(/;)2.(l + 0(p)) 

^ (^ln2 + 0(n- 1)) . (A.15) 

Finally we get 



n —\ 

because (see Eqs. f|TI)|) and (fT^ ) 



1 +0(n-l))/2(n), (A.16) 



^:||^ ~ Jo-V) « • (A.17) 
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Thus, the result of ()A.16|) does not depend on the specific value of A. Anal- 
ogously we find 

"1/^ - dy 



e/2 

M ( 1 



/ -/,%) = / ^/|(y) + ln2/|(x)-(l + 0(p)^ 

Jxi2 y Jx y ^ ^ 

+ ln2 + 0(n-l) ) /2H • (A.18) 



n — 1 

By comparing Eq. ()A.18|) for A=2 and Eq. pij). we get 

2"-i /2(n|l/2) = + ln2 + 0(n - 1)) Mn) (A.19) 

n — 1 \n — 1 y 

and conclude that 

Mn\l/2) = Mn) + 0{in-l)'). (A.20) 

This completes the proof of Eq. (j2Z|). Equation ()A.20|) is very important and 
helps us to sum the regular parts of Fag{x). We should remember, however, 
that this result holds only for x-independent input distributions, that is for 
Bg = Bg = Om Eq. ©. 

Comparing Eqs. ()A.4|) . ()A.11|) and ()A.15|) . we get finally 

r ^/|(y) = ln2/|(x)-(l + 0(p)) ^ ln2/2H + 0(n-l). (A.21) 

Jx/2 y ^ ^ 

B Appendix: The coefficients C^^^ 

Here we give the explicit values of the coefficients Cli^{n = 1) = C*^"^ ap- 
pearing in Eqs. (jH^ and (j^ : 

,__ 1813 297 

"5 ~ 2880 ' "5 ~ 2880 32 ^ ' 

297 



_ 283 ^2 — 33 1813 
~ ~6480 ' 99 - ~ 



99 32 \^ 4374 y ^ ) 
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For the coefficient C'^'^ one must consider separately the singular and 
regular part, 



Ci+ = + = 1) 



aa 



n — 1 



where 



'-'99 



2673 
~32~ 



516577 
5184 



(B.2) 



C Appendix: Proof of the simplified form for 
the functions Fi and F2 

The validity of the simplified form for the functions Fi and F2, given in 
Eqs. (|56 p -(|6i p . is based on the following estimates. As in Appendix A, we 
consider separately the three components of the functions Fi and F2 given 
in Eq. (IKHll . 

a) Consider first the " " component. The function (/~(x,M^))^ is x- 

independent since we assumed Bg = Bg = : (/^"(x, M^))^ = {fg{M'^)Y, 
so 



frin,M') 



n - 1 

Then (see Eqs. and (0) we have 



(C.l) 



'a 

M 



n — 1 



2nx2 M d+v/{n-l) 



(n-1) 



(C.2) 



2 ^■'9 



M- 



-1 1 



;ji(a)(/-(x,M2))2 



Thus, 



{A-fUa)-e-''-^'^-(l + 0{p)) , 
(A^-)2iji(a)-e-2'^-«-(l + 0(p)) . 



(C.3) 

(C.4) 
(C.5) 
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b) Consider now the "H — " component. The function /^(x, M^) can be 
represented as 



Note that 

J n - 

1 t . . , J r-l 

Io{a) . (C.7) 



/q(o-) ^ ^d+v/(n-l) 1 ^-d+w/(n-l) 

n — 1 



n — 1 

In the same way we obtain 



Ma) ® Jo(a) ^ ^e<l+^/(n-i)^g-dW(n-i) 

n — 1 n — 1 

1 ^-dV/(n-i) i^l;^ (C.8) 



{n — 1)2 p 

Thus, we get 

F+- = -pJi(a)® (/;(x,M2)/;(x,M 

= A^-A+/o(a)-e-('^-«+^+«)'"(l + 0(p)) , (C.9) 
F/- = Jo(a)® (/;(x,M2)/-(x,M2)) 

= A;v4+^/i(a) ■e-('^-«+'^+«)'"(^l + 0(p2)j . (C.IO) 

c) Finally, let us consider the "++" component. The problem of finding the 
small-x behavior of Jo (5") Io{cr) can be solved as follows. From the tables 
of the inverse Mellin transforms in Ref. EHl, we obtain 



Jo 

where y = — log a; and Re z/ > —1, Re 7 > 0. 
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Equation ()C.11|) gives, for /? = z/ = 0, a = —MggW and 7 = n — 1, the 
Mellin transform of /gl*^)- 

n — 1 \ n — 1 I 

The same result can be obtained in an alternative way by performing the 
Mellin transform of the series 

, 1 ^ (a)^-r(m + l/2) 

From this first step we get 

M^)®im ^ -^e-^V(-^*)/(-^)/of^^) (C.13) 
(?7, — I)"' \ n — 1 / 

and the inverse Mellin transform of the right hand side of ()C.13|) will be the 
final answer for F2~^ (apart from the overall factor (A^)"^). 

We consider now Eq. (frrnl for z/ = 0. By taking the sum of the derivative 
of Eq. ()C.11|) with respect to a with its derivative with respect to (3, we obtain 
the important formula 

1 



7^ \ 27 

where, as before, 7 = n — 1 and y = — logx. With the substitutions 

a ^ -d+{3w + t), /3^-d+{t-w), (C.15) 

Eq. ()C.14|) provides the final answer, since the coefficient of x"^"^ in the 
integral at the left hand side gives the required inverse Mellin transform. 
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A similar procedure can be applied to determine pJiip-) ® Io{cr), which 
is needed to calculate F^^, starting from the formula ()C.11|) . with the same 
values of the parameters a and (3 as in ()C.15|) and y = — logx. Notice that 
P < and hence Z?^/^ is pure imaginary. 

After introducing, for the sake of simplicity, the notation z 
i \(3^^'^\, we get 



a 



1/2 



+ 



6{l-x)-pJi{a)\ ® (^/+(x,M2) 
(A+)2/o(V/')/o(V/')-e-2^+«"'-fl + 0{p)), (C.16) 



while, with the same definition of the variables a, (3 and we get 



y'" [zh{zy'l')h{zy'") 



- ^/o(V/')/i(zV/^)l -e-^^+W-- fl + 0{p')) . (C.17) 



This completes the calculation of the simplified form of the functions Fi 
and F2. 
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